New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects.

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

Sexual stage-specific A-to-I mRNA editing is mediated by tRNA-editing enzymes in fungi.

A-to-I RNA editing catalyzed by adenosine-deaminase-acting-on-RNA (ADARs) was assumed to be unique to metazoans because fungi and plants lack ADAR homologs. However, genome-wide messenger RNA (mRNA) editing was found to occur specifically during sexual reproduction in filamentous ascomycetes. Because systematic characterization of adenosine/cytosine deaminase genes has implicated the involvement of TAD2 and TAD3 orthologs in A-to-I editing, in this study, we used genetic and biochemical approaches to characterize the role of FgTAD2, an essential adenosine-deaminase-acting-on-tRNA (ADAT) gene, in mRNA editing in Fusarium graminearum. FgTAD2 had a sexual-stage-specific isoform and formed heterodimers with enzymatically inactive FgTAD3. Using a repeat-induced point (RIP) mutation approach, we identified 17 mutations in FgTAD2 that affected mRNA editing during sexual reproduction but had no effect on transfer RNA (tRNA) editing and vegetative growth. The functional importance of the H352Y and Q375*(nonsense) mutations in sexual reproduction and mRNA editing were confirmed by introducing specific point mutations into the endogenous FgTAD2 allele in the wild type. An in vitro assay was developed to show that FgTad2-His proteins purified from perithecia, but not from vegetative hyphae, had mRNA editing activities. Moreover, the H352Y mutation affected the enzymatic activity of FgTad2 to edit mRNA but had no effect on its ADAT activity. We also identified proteins co-purified with FgTad2-His by mass spectrometry analysis and found that two of them have the RNA recognition motif. Taken together, genetic and biochemical data from this study demonstrated that FgTad2, an ADAT, catalyzes A-to-I mRNA editing with the stage-specific isoform and cofactors during sexual reproduction in fungi.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants.

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

Programmable RNA base editing with photoactivatable CRISPR-Cas13.

CRISPR-Cas13 is widely used for programmable RNA interference, imaging, and editing. In this study, we develop a light-inducible Cas13 system called paCas13 by fusing Magnet with fragment pairs. The most effective split site, N351/C350, was identified and found to exhibit a low background and high inducibility. We observed significant light-induced perturbation of endogenous transcripts by paCas13. We further present a light-inducible base-editing system, herein called the padCas13 editor, by fusing ADAR2 to catalytically inactive paCas13 fragments. The padCas13 editor enabled reversible RNA editing under light and was effective in editing A-to-I and C-to-U RNA bases, targeting disease-relevant transcripts, and fine-tuning endogenous transcripts in mammalian cells in vitro. The padCas13 editor was also used to adjust post-translational modifications and demonstrated the ability to activate target transcripts in a mouse model in vivo. We therefore present a light-inducible RNA-modulating technique based on CRISPR-Cas13 that enables target RNAs to be diversely manipulated in vitro and in vivo, including through RNA degradation and base editing. The approach using the paCas13 system can be broadly applicable to manipulating RNA in various disease states and physiological processes, offering potential additional avenues for research and therapeutic development.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns.

BACKGROUND: Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events. RESULTS: In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species. CONCLUSIONS: Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Narrowing down the candidates of beneficial A-to-I RNA editing by comparing the recoding sites with uneditable counterparts.

Adar-mediated adenosine-to-inosine (A-to-I) RNA editing mainly occurs in nucleus and diversifies the transcriptome in a flexible manner. It has been a challenging task to identify beneficial editing sites from the sea of total editing events. The functional Ser>Gly auto-recoding site in insect Adar gene has uneditable Ser codons in ancestral nodes, indicating the selective advantage to having an editable status. Here, we extended this case study to more metazoan species, and also looked for all Drosophila recoding events with potential uneditable synonymous codons. Interestingly, in D. melanogaster, the abundant nonsynonymous editing is enriched in the codons that have uneditable counterparts, but the Adar Ser>Gly case suggests that the editable orthologous codons in other species are not necessarily edited. The use of editable versus ancestral uneditable codon is a smart way to infer the selective advantage of RNA editing, and priority might be given to these editing sites for functional studies due to the feasibility to construct an uneditable allele. Our study proposes an idea to narrow down the candidates of beneficial recoding sites. Meanwhile, we stress that the matched transcriptomes are needed to verify the conservation of editing events during evolution.

ADAR-mediated RNA editing regulates PVR immune checkpoint in colorectal cancer.

Recent studies have revealed that tumor immunotherapy resistance is influenced by ADAR-mediated RNA editing, but its targets remain unelucidated. Our current study identified the poliovirus receptor (PVR) oncogene, which encodes an immune checkpoint in colorectal cancer (CRC), as a potential target for RNA editing. We performed transcriptome sequencing analysis and experimental validation in two Chinese CRC cohorts. PVR and ADAR expressions significantly increased in CRC tumors and showed positive correlations in both cohorts, coupled with upregulated PVR RNA editing in CRC tumors. Manipulation of ADAR expression by over-expression or knockdown substantially changed PVR expression and RNA editing in HTC116 CRC cells. Luciferase reporter and actinomycin D assays further revealed that RNA editing in PVR 3'-UTR could upregulate PVR RNA expression, probably by increasing the RNA stability. By increasing PVR expression, ADAR-mediate RNA editing might contribute to tumor- and immune-related gene functions and pathways in CRC. Moreover, a signature combining PVR RNA editing and expression showed promising predictive performance in CRC diagnosis in both Chinese CRC cohorts. Our findings thus highlight the importance of ADAR-mediated RNA editing in PVR up-regulation in CRC tumors and provide new insight into the application of PVR RNA editing as a novel diagnostic biomarker for CRC.

Adaptive advantages of restorative RNA editing in fungi for resolving survival-reproduction trade-offs.

RNA editing in various organisms commonly restores RNA sequences to their ancestral state, but its adaptive advantages are debated. In fungi, restorative editing corrects premature stop codons in pseudogenes specifically during sexual reproduction. We characterized 71 pseudogenes and their restorative editing in Fusarium graminearum, demonstrating that restorative editing of 16 pseudogenes is crucial for germ tissue development in fruiting bodies. Our results also revealed that the emergence of premature stop codons is facilitated by restorative editing and that premature stop codons corrected by restorative editing are selectively favored over ancestral amino acid codons. Furthermore, we found that ancestral versions of pseudogenes have antagonistic effects on reproduction and survival. Restorative editing eliminates the survival costs of reproduction caused by antagonistic pleiotropy and provides a selective advantage in fungi. Our findings highlight the importance of restorative editing in the evolution of fungal complex multicellularity and provide empirical evidence that restorative editing serves as an adaptive mechanism enabling the resolution of genetic trade-offs.

REDH: A database of RNA editome in hematopoietic differentiation and malignancy.

BACKGROUND: The conversion of adenosine (A) to inosine (I) through deamination is the prevailing form of RNA editing, impacting numerous nuclear and cytoplasmic transcripts across various eukaryotic species. Millions of high-confidence RNA editing sites have been identified and integrated into various RNA databases, providing a convenient platform for the rapid identification of key drivers of cancer and potential therapeutic targets. However, the available database for integration of RNA editing in hematopoietic cells and hematopoietic malignancies is still lacking. METHODS: We downloaded RNA sequencing (RNA-seq) data of 29 leukemia patients and 19 healthy donors from National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, and RNA-seq data of 12 mouse hematopoietic cell populations obtained from our previous research were also used. We performed sequence alignment, identified RNA editing sites, and obtained characteristic editing sites related to normal hematopoietic development and abnormal editing sites associated with hematologic diseases. RESULTS: We established a new database, "REDH", represents RNA editome in hematopoietic differentiation and malignancy. REDH is a curated database of associations between RNA editome and hematopoiesis. REDH integrates 30,796 editing sites from 12 murine adult hematopoietic cell populations and systematically characterizes more than 400,000 edited events in malignant hematopoietic samples from 48 cohorts (human). Through the Differentiation, Disease, Enrichment, and knowledge modules, each A-to-I editing site is systematically integrated, including its distribution throughout the genome, its clinical information (human sample), and functional editing sites under physiological and pathological conditions. Furthermore, REDH compares the similarities and differences of editing sites between different hematologic malignancies and healthy control. CONCLUSIONS: REDH is accessible at http://www.redhdatabase.com/ . This user-friendly database would aid in understanding the mechanisms of RNA editing in hematopoietic differentiation and malignancies. It provides a set of data related to the maintenance of hematopoietic homeostasis and identifying potential therapeutic targets in malignancies.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target.

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Improvement of C-to-U RNA editing using an artificial MS2-APOBEC system.

RNA cytidine deamination (C-to-U editing) has been achieved using the MS2-apolipoprotein B-editing catalytic polypeptide-like (APOBEC)1 editing system. Here, we fused the cytidine deaminase (CDA) enzymes APOBEC3A and APOBEC3G with the MS2 system and examined their RNA editing efficiencies in transfected HEK 293T cells. Given the single-stranded RNA preferences of APOBEC3A and APOBEC3G, we designed unconventional guide RNAs that induced a loop at the target sequence, allowing the target to form a single-stranded structure. Because APOBEC3A and APOBEC3G have different base preferences (5'-TC and 5'-CC, respectively), we introduced the D317W mutation into APOBEC3G to convert its base preference to that of APOBEC3A. Upon co-transfection with a guide RNA that induced the formation of a 14 nt loop on the target sequence, MS2-fused APOBEC3A and APOBEC3G showed high editing efficiency. While the D317W mutation of APOBEC3G led to a slight improvement in editing efficiency, the difference was not statistically significant. These findings indicate that APOBEC3A and APOBEC3G can induce C-to-U RNA editing when transfected with a loop guide RNA. Moreover, the editing efficiency of APOBEC3G can be enhanced by site-specific mutation to alter the base preference. Overall, our results demonstrate that the MS2 system can fuse and catalyze reactions with different enzymes, suggesting that it holds an even greater potential for RNA editing than is utilized currently.

Profiling of RNA editing events in plant organellar transcriptomes with high-throughput sequencing.

RNA editing is a crucial post-transcriptional modification process in plant organellar RNA metabolism. rRNA removal-based total RNA-seq is one of the most common methods to study this event. However, the lack of commercial kits to remove rRNAs limits the usage of this method, especially for non-model plant species. DSN-seq is a transcriptome sequencing method utilizing duplex-specific nuclease (DSN) to degrade highly abundant cDNA species especially those from rRNAs while keeping the robustness of transcript levels of the majority of other mRNAs, and has not been applied to study RNA editing in plants before. In this study, we evaluated the capability of DSN-seq to reduce rRNA content and profile organellar RNA editing events in plants, as well we used commercial Ribo-off-seq and standard mRNA-seq as comparisons. Our results demonstrated that DSN-seq efficiently reduced rRNA content and enriched organellar transcriptomes in rice. With high sensitivity to RNA editing events, DSN-seq and Ribo-off-seq provided a more complete and accurate RNA editing profile of rice, which was further validated by Sanger sequencing. Furthermore, DSN-seq also demonstrated efficient organellar transcriptome enrichment and high sensitivity for profiling RNA editing events in Arabidopsis thaliana. Our study highlights the capability of rRNA removal-based total RNA-seq for profiling RNA editing events in plant organellar transcriptomes and also suggests DSN-seq as a widely accessible RNA editing profiling method for various plant species.

Genome-Wide Analysis on Driver and Passenger RNA Editing Sites Suggests an Underestimation of Adaptive Signals in Insects.

Adenosine-to-inosine (A-to-I) RNA editing leads to a similar effect to A-to-G mutations. RNA editing provides a temporo-spatial flexibility for organisms. Nonsynonymous (Nonsyn) RNA editing in insects is over-represented compared with synonymous (Syn) editing, suggesting adaptive signals of positive selection on Nonsyn editing during evolution. We utilized the brain RNA editome of Drosophila melanogaster to systematically study the LD (r2) between editing sites and infer its impact on the adaptive signals of RNA editing. Pairs of editing sites (PESs) were identified from the transcriptome. For CDS PESs of two consecutive editing sites, their occurrence was significantly biased to type-3 PES (Syn-Nonsyn). The haplotype frequency of type-3 PES exhibited a significantly higher abundance of AG than GA, indicating that the rear Nonsyn site is the driver that promotes the editing of the front Syn site (passenger). The exclusion of passenger Syn sites dramatically amplifies the adaptive signal of Nonsyn RNA editing. Our study for the first time quantitatively demonstrates that the linkage between RNA editing events comes from hitchhiking effects and leads to the underestimation of adaptive signals for Nonsyn editing. Our work provides novel insights for studying the evolutionary significance of RNA editing events.

Extracellular vesicle­mediated RNA editing may underlie the heterogeneity and spread of hepatocellular carcinoma in human tissue and in vitro.

Extracellular vesicles (EVs) produced by various cells, including tumor cells, carry biomolecules to neighboring cells. In hepatocellular carcinoma (HCC), adenosine to inosine RNA editing of antizyme inhibitor 1 (AZIN1), specifically regulated by adenosine deaminase acting on RNA­1 (ADAR1), promotes carcinogenesis. The present study examined if EVs and ADAR1 in the EVs released from HCC cells are transferred to neighboring cells in co­culture systems and reporter assay. Distribution of the ADAR1 expression in human tissues were examined by immunohistochemistry. EVs released from HCC cells containing ADAR1 were delivered to neighboring HCC cells and non­cancerous hepatocytes. The increased ADAR1 protein levels resulted in serine to glycine substitution at residue 367 of AZIN1, which augmented transformation potential and increased aggressive behavior of cancer cells. In clinically resected samples, ADAR1 distribution was highly heterogeneous within the tumor specimen and denser in non­cancerous tissue surrounding the HCC tissue. These observations suggested that ADAR1 protein may be delivered from HCC cells to neighboring cells via EVs and that EV­mediated RNA editing may serve a pivotal role in determining HCC heterogeneity and spread.

A Genome-Wide Analysis of the Pentatricopeptide Repeat Protein Gene Family in Two Kiwifruit Species with an Emphasis on the Role of RNA Editing in Pathogen Stress.

Kiwifruit is a perennial fruit tree with high nutritional and economic value; however, various pathogen stresses have resulted in reductions in its yield and quality. Pentatricopeptide repeat proteins (PPRs), characterized by tandem repetitions of 35 amino acid motifs, play roles in RNA editing, mRNA stability, and splicing. They may also regulate plant development and growth. Nevertheless, the roles of PPRs in plant development and disease resistance remain unclear. In this study, we focused on the roles of PPRs in the fruit development and pathogen stress of kiwifruit and conducted a series of analyses of the PPR gene family in two representative kiwifruit species (Actinidia chinensis (Ach) and Actinidia eriantha (Ace)) with markedly different degrees of disease resistance. A total of 497 and 499 PPRs were identified in Ach and Ace, respectively. All the kiwifruit PPRs could be phylogenetically divided into four subfamilies. There were about 40.68% PPRs predicted to be localized to mitochondria or chloroplasts. A synteny analysis showed that the expansion of the kiwifruit PPRs mainly originated from segmental duplication. Based on RNA-seq data from the fruit over 12 periods of development and maturity, a weighted correlation network analysis suggested that two PPRs, Actinidia20495.t1 and Actinidia15159.t1, may be involved in fruit development and maturation. In addition, we observed different responses with respect to the expression of PPRs and RNA editing between resistant and susceptible kiwifruits following infection with pathogenic bacteria, indicating the regulatory role of PPRs in the stress response via the modulation of RNA editing. The differentially expressed upstream transcription factors of the PPRs were further identified; they may regulate resistance adaption by modulating the expression of the PPRs. Collectively, these results suggest that PPRs play roles in the development and disease resistance of kiwifruit and provide candidate genes for further clarifying the resistance mechanisms in kiwifruits.

Examination of Factors Affecting Site-Directed RNA Editing by the MS2-ADAR1 Deaminase System.

Adenosine deaminases acting on RNA (ADARs) have double-stranded RNA binding domains and a deaminase domain (DD). We used the MS2 system and specific guide RNAs to direct ADAR1-DD to target adenosines in the mRNA encoding-enhanced green fluorescence protein. Using this system in transfected HEK-293 cells, we evaluated the effects of changing the length and position of the guide RNA on the efficiency of conversion of amber (TAG) and ochre (TAA) stop codons to tryptophan (TGG) in the target. Guide RNAs of 19, 21 and 23 nt were positioned upstream and downstream of the MS2-RNA, providing a total of six guide RNAs. The upstream guide RNAs were more functionally effective than the downstream guide RNAs, with the following hierarchy of efficiency: 21 nt > 23 nt > 19 nt. The highest editing efficiency was 16.6%. Off-target editing was not detected in the guide RNA complementary region but was detected 50 nt downstream of the target. The editing efficiency was proportional to the amount of transfected deaminase but inversely proportional to the amount of the transfected guide RNA. Our results suggest that specific RNA editing requires precise optimization of the ratio of enzyme, guide RNA, and target RNA.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species.

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

A-to-I RNA editing shows dramatic up-regulation in osteosarcoma and broadly regulates tumor-related genes by altering microRNA target regions.

A-to-I RNA editing is a prevalent type of RNA modification in animals. The dysregulation of RNA editing has led to multiple human cancers. However, the role of RNA editing has never been studied in osteosarcoma, a complex bone cancer with unknown molecular basis. We retrieved the RNA-sequencing data from 24 primary osteosarcoma patients and 3 healthy controls. We systematically profiled the RNA editomes in these samples and quantitatively identified reliable differential editing sites (DES) between osteosarcoma and normal samples. RNA editing efficiency is dramatically increased in osteosarcoma, presumably due to the significant up-regulation of editing enzymes ADAR1 and ADAR2. Up-regulated DES in osteosarcoma are enriched in 3'UTRs. Strikingly, such 3'UTR sites are further enriched in microRNA binding regions of gene EMP2 and other oncogenes, abolishing the microRNA suppression on target genes. Accordingly, the expression of these tumor-promoting genes is elevated in osteosarcoma. There might be an RNA editing-dependent pathway leading to osteosarcoma. We expanded our knowledge on the potential roles of RNA editing in oncogenesis. Based on these molecular features, our work is valuable for future prognosis and diagnosis of osteosarcoma.

DYW cytidine deaminase domains have a long-range impact on RNA recognition by the PPR array of chimeric plant C-to-U RNA editing factors and strongly affect target selection.

The protein factors for the specific C-to-U RNA editing events in plant mitochondria and chloroplasts possess unique arrays of RNA-binding pentatricopeptide repeats (PPRs) linked to carboxy-terminal cytidine deaminase DYW domains via the extension motifs E1 and E2. The E1 and E2 motifs have distant similarities to tetratricopeptide repeats known to mediate protein-protein interactions but their precise function is unclear. Here, we investigate the tolerance of PPR56 and PPR65, two functionally characterized RNA editing factors of the moss Physcomitrium patens, for the creation of chimeras by variably replacing their C-terminal protein regions. Making use of a heterologous RNA editing assay system in Escherichia coli we find that heterologous DYW domains can strongly restrict or widen the spectrum of off-targets in the bacterial transcriptome for PPR56. Surprisingly, our data suggest that these changes are not only caused by the preference of a given heterologous DYW domain for the immediate sequence environment of the cytidine to be edited but also by a long-range impact on the nucleotide selectivity of the upstream PPRs.